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ABSTRACT 



Aims. We present a new semi-analytical model of galaxy formation, GECO (Galaxy Evolution COde), aimed at a better understanding 
of when and how the two processes of star formation and galaxy assembly have taken place, by comparison with a wide variety of 
recent data on the evolutionary galaxy mass functions and star-formation histories. 

Methods. Our model is structured into a Monte Carlo algorithm based on the Extended Press-Schechter theory, for the representation 
of the merging hierarchy of dark matter halos, and a set of analytic algorithms for the treatment of the baryonic physics, including 
classical recipes for the gas cooling, the star formation time-scales, galaxy mergers and SN feedback. Together with the galaxies, the 
parallel growth of BHs is followed in time and their feedback on the hosting galaxies is modelled. We set the model free parameters 
by matching with data on local stellar mass functions and the BH-bulge relation at z = 0. 

Results. Based on such local boundary conditions, we investigate how data on the high-redshift universe constrain our understanding 
of the physical processes driving the evolution, focusing in particular on the assembly of stellar mass and on the star formation history. 
Since both processes are currently strongly constrained by cosmological near- and far-IR surveys with the Spitzer Space Telescope, 
the basic physics of the ACDM hierarchical clustering concept of galaxy formation can be effectively tested by us by comparison 
with the most reliable set of observables using a minimal number of free parameters. 

Conclusions. Our investigation shows that when the time-scales of the stellar formation and mass assembly are studied as a function 
of dark matter halo mass and the single galaxy stellar mass, the 'downsizing' fashion of star formation appears to be a natural outcome 
of the model, reproduced even in the absence of the AGN feedback. On the contrary, the stellar mass assembly history turns out to 
follow a more standard hierarchical pattern progressive in cosmic time, with the more massive systems assembled at late times mainly 
through dissipationless mergers. 

Key words, galaxies: evolution - galaxies: formation - galaxies: halos. 



1. Introduction 

In the past decade, several observational evidences were accu- 
mulating in favour of the ACDM paradigm for structure for- 
mation, now quite a successful rendition of the hierarchical 
clustering s cenario for cosmic struc ture formation. In its stan- 
dard form jBlumenthal et al.L [1984'), it predicts that structures 
formed from primordial fluctuations of the density field am- 
plified during inflation and then collapsed to form the virial- 
ized structures that we see nowadays. The most compelling 
support in favour of this paradigm comes from the mea- 
surements of the temperature anisotropics of the cosmic mi- 
crowave background (Spergel et al., 2003, 2007). Further ev- 
idences are due to the measurements of the power spectrum 
of galaxy clu stering f rom large surveys o f the local universe 
dPercival et al ., 2002; T egmark et al.L l2004h . the evidence for an 
accelerated expansion of the universe as i nferred from high- 
redshift type la supernovae obse rvations dRiess et al.L Il998t 
iPerlmutter. Turner. & White , 1999) and th e baryon fraction ob- 
served in rich clusters dWhite et al Lll993h . 

Devised to study galaxy evolution within this cosmologi- 
cal framework, the semi-analytical approach favours a relatively 
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simple handling of the main physical parameters and under- 
standing of their possible role in drivin g the evolution. Thi s 
modelling has its root in the work of White & ReesI ( Il978l) . 
where it was proposed that galaxy formation is a two-stage pro- 
cess, with dark matter halos forming in a dissipationless grav- 
itational collapse and galaxies forming inside them following 
the radiativ e cooling of baryons. A lthough White & ReesI (1 19781) 
and, after, IWhite & FrenkI d 199 ID . based their work only o n 
the analytic Press-Schechter theory dPress & SchechteR 1 19741) . 
predicting only average quantities, subsequently a number of 
works followed their prescriptions using Monte Carlo (MC) 
merger trees. The MC approach allows to obtain several real- 
izations of the merging history of individual dark mat ter ha - 
los. This approach was pi o neered by [Lacey &Coig (fT99l : 
iKauffmann & White! dl993l) : ICole e^ al.' (1994 ), and th e n fol- 
lowed by a number of authors ((Somerville & Kolatl 119991 : 
IShefli & Lernsoiil 119991; IZentneii |2007|). The great advantage of 
the semi-analytical method (SAM), apart from being computa- 
tionally very fast and flexible, is the fact that it is possible to 
compute merging histories with arbitrary mass resolution. The 
alternative approach that can be followed in order to track the 
evolution of the dark matter halos, is through the use of large 
cosmological N-body simulations. Their large computational re- 
quirements are compensated by the amount of information that 
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can be achieved. For instance, in the Millennium Simulation 
(ISpringel. Yoshidar& Whitel 1200 Ih . the evolution of substruc- 
tures in massive halos can be followed in time, with the re- 
sults of a more detailed information about the galaxy dynam- 
ics a nd the influence of th e cosmic environment on the pro- 
cess (iDe Lucia et al.L |2004|) . In the literature, various examples 
of this "hybrid" approach, which make use of N-body simula- 
tion for the dark matter evolution and the SAM technique for 
the baryonic physics, have been publishe d dHatton et al.L 120031: 
iDe Lucia et a"Lll2006HCro'ton et al.Ll2006l) . 

In the present work we employ a MC merger tree, mainly 
because this allowed us to test the parameter space of the semi- 
analytical model with much more flexibility than using N-body 
simulations, and allows to compute merging histories down to 
arbitrary low mass resolution. 

While the treatment of the evolution of dark mat- 
ter structures is relatively simple, as being determined 
only by gravity, the physics of the baryons inside ha- 
los is much more com plex to describe. In the most 

classical models ('KaufFmann. White. & Guiderdoni', 

Baugh. Cole. & Frenk, 199 6; Kauffma nn et aL , 
omerville & Primackj 1 1999b ICole et all l200ot 
Menci & Cavaliere^ ,2002) . gas cooling, star formation, SN 
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feedback and galaxy mergers are included. In recent years, it has 
become clear that some other form of highly energetic feedback 
is needed to prevent star formation in massive galaxies at recent 
epochs, where the SN feedback is ineffective. The need of such 
form of feedback is required in order to avoid the overcooling in 
massive halos, and hence the overabundance of galaxies at the 
bright-end of the luminosity function. This source of feedback 
is commonly found in the AGN energy production. This effect, 
supported by the observational findings o f a tight correlation 
between the BH and the bulge size (iFerrarese & Merriti , 
I2OOOI; paring &Rix, 2004^ was imple mented in different 
way s by Kauffmann & Haehnelt (2000); Bower etal. (2006^); 
ICroton et al. (2006) ; Menci et al.. (2006); Somerville et alj 
(I2008h . 



Alternative to the AGN feedback, the shutdown of star for- 
mation above a critical halo mass has been implemented as 
a quenchin g mechanism in massive galaxies (ICattaneo et al.L 
I2OO6II2OO8I) . motivated by the prediction of stable shock heating 
for ha los more massive than this threshold (iDekel & Birnboiml 
I2006h . 

Following the prescriptions of these models, we have built 
a new semi-analytical models, the Galaxy Evolution COde 
(GECO), whose aim is to identify a few key physical parame- 
ters and modify them by comparing with several basic properties 
of the galaxy population at z = 0, as well as at high redshifts. 
Our main observational reference in this paper is the redshift- 
dependent stellar mass function of galaxies, which, from a suit- 
able choice for the stellar Initial Mass Function (IMF), is a robust 
descriptor of the star formation history and the mass assembly 
history of galaxies. 

The structure of the paper is as follows. In ^we describe 
the Monte Carlo merger tree used in the model. In ^we in- 
troduce the ingredients of the baryonic model. In ^we explain 
how the free parameters are set and provide a table for them. In 
^ the basic results for the local universe are presented, while 
in ^ we focus on the high redshift predictions. We conclude 
in I?] Throughout the rest of the paper we assume a "concor- 
dance" cosmological model, with Q,„ - 0.3, Q.a = 0.7, h - 0.7, 
erg - 0.9 and n - \ (power index of the primordial power 
spectrum). However, when needed, we also show the depen- 
dence on the cosmological parameters showing the results for 



the new WMAP5 dataset, which are: Q.^ = 0.258, = 742, 
h = 0.719, crs = 0.796 and n = 0.963 (iDunklev et al.Ll2009h . 

2. Merger Tree 

2.1. The formalism 

In the hierarchical theory of structure formation the density con- 
trast 6 - (p - Pb)/Pb, where pi, is the mean density of the uni- 
verse and p is the perturbation density, is a Gaussian random 
field that at early times grows linearly. When the density con- 
trast reaches a critical value 6c, given by the spherical collapse 
theory, the perturbation starts to collapse and form virialized ha- 
los. In order to associate masses to these collapsed objects, the 
Press-Schechter approach considers the fluctuations whose den- 
sity contrast smoothed on a scale Rm, through a spatial window 
function W{r,RM), exceeds 6^. 

The Gaussian random variable S has zero mean and variance 
S{M) = cP-(M) which is linked to the mass through the Power- 
spectrum P{k): 



P(k)W^(kr)k^dk, 



(1) 



where W^{kr) is the Fourier transform of the spatial win- 
dow function. In the excursion-set approach, developed by 
iBond et al ] (1199 Ih . the value of 6 executes a random walk as 
the smoothing scale Rm (or M) is changed. We can consider the 
trajectories in the plane (S, 6c) and associate the fraction of mat- 
ter in collapsed objects in the range dM around M at the time 
t with the fraction of trajectories that make the first upcrossing 
through the threshold 6c in the interval S,S + dS . This results in 
the well-known PS mass function: 



fiS,6c)dS = 



1 



Sc 



■ exp 



(-1)- 



(2) 



(iPress & Schechteiill974 . 

In the Extended Press-Schechter model (iBond et al.L 119911 
iLacev & Cole, 1993) it is also possible to derive the conditional 
mass function, the fraction of trajectories in halos with mass M\ 
at z\ that are in halos with mass Mq > Mi at zo < Zi : 



f{SuSci\So,6co)dSi = 
1 (6ci - 6co) 



V2^(5i-5„)3/2 



-exp 



2{Si-So) 



dS 1 



(3) 



Converting from mass weighting to number weighting we 
obtain the average number of progenitors at z\ in the mass inter- 
val dM[ around Mi which at redshift zo has merged to form an 
Mq halo: 



dP 
Mo 



(Mi,zi\Mo,Zo)dMi = 
dS 



— f{S 1 , 6ci \Sq, 6co)\ IdMi 

Ml \dM\ 



(4) 



2.2. The partition algorithm 

For generating Monte Carlo realizations of the merging history 
of dark matter halos we u se the partition algorithm described 
bv'Sheth & LemsonI (Il999h . hereafter SL99, which is exact for 
a white-noise power-spectrum, but needs some modifications in 
order to be applied to a ACDM fluctuation spectrum. Since we 
will describe in detail the algorithm and the test done to probe 
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its consistency in a separate paper, we give here only a brief 
description of it. 

The algorithm is based on the assumption that for a white- 
noise spectrum mutually disconnected volumes are mutually in- 
dependent, as it is analytically demonstrated in the Appendix of 
SL99. Let us call parent the halo of mass Mo which exists at z-Q 
and progenitors the smaller halos that derive from it at higher 
redshift. Suppose to partition an halo of mass Mq into progeni- 
tors by choosing first one progenitor and then an other one from 
the remaining mass and so on until the remaining mass falls be- 
low a certain minimum mass that is our mass resolution Mres- We 
have a probability of finding the first progenitor of mass Mi , cor- 
responding to 5 1, given by f{S\,6c\\So,6ca) (equation[3]) and we 
choose a mass drawing a random number from this distribution. 
We can consider the halo Mq as a region of size Vm^ - Mobo/pi,, 
where /jQ = 1/(1 -H(Jo)- If the first progenitor has mass Mi it occu- 
pies a volume Vm, = Mi/ji/p/,. Theremainingmass/? = M^-Mi 
is distributed in the volume Vr - RbR/pt,. Hence, for the conser- 
vation of volumes, the density in this region is given by 



3eR 



del 



I -{Ml /Mo) 



(5) 



to lowest order in 6 (see equation 5 of SL99). Now, the second 
progenitor must be chosen from the remaining mass R and the 
probability that it has mass M2, corresponding to S2, is given 
by f(S 2, Sci\S (R), 6cr)- Therefore we choose a random number 
from this distribution. Then we iterate the process by continuing 
to find progenitors until the remaining mass is below the mass 
resolution. 

The assumption that disconnected volumes are mutually in- 
dependent is right only for white noise power spectrum, but we 
want to build merger trees for a ACDM spectrum. SL99 show 
that by applying the same algorithm for scale-free spectra with 
n 5t leads to inconsistencies in the excursion set mean val- 
ues. Since we know that all the excursion set quantities are in- 
dependent of the power spectrum when they are expressed in 
terms of the variance rather than the mass, we run the algo- 
rithm for the white-noise case. Then we consider each chosen 
M,„, not as a progenitor with mass M,„,, but as a region pop- 
ulated by some number v of objects having mass M^cdm- In 
practise, we normalize the white-noise power spectrum in such 
a way that the mass variance of the two spectra is the same in 
Mq: S[) = S ACDMiMo). Then, we choose Si from the distribu- 
tion of equation [3] We associate to this variance a mass given 
by the white-noise spectrum M,„, - SoMq/S i and one given by 
the ACDM spectrum M/^cdm- We consider a number v of ha- 
los with this mass, given by v = Mt„,/MACDM, approximated to 
the nearest integer. Then, we iterate the procedure on each new 
progenitor 

As a test of the consistency of our MC merger tree with the 
PS predictions, in Figure [T] we show the global mass function 
of dark matter halos. To compute it, we use a grid of 35 par- 
ent halo masses, ranging from lO'^M© to lO'^M©, logarithmi- 
cally spaced, and simulate 100 realizations for each parent halo. 
Then we weight each halo with the PS number density at z=0. 
The results do not depend much on the number of parent halo 
masses used, but they are quite sensitive to the number of real- 
izations for each halo. The 100 realizations used turn out to be 
a good compromise between the good sampling desired and the 
computational-time required. We build a grid in time of 52 inter- 
vals logarithmically spaced in redshift, ranging from z=20 when 
the time-step is approximately equal to 0.05 Gyr to z=0 with 
a ~0.8 Gyr time-step. Because our method for building merger 




log(M) 

Fig. 1. Global mass function for dark matter halos at different 
redshifts. Histograms represent the mass function derived from 
our MC merger tree, while solid lines show the PS predictions. 



tree is insensitive to the time-step used, the choice of the time- 
step grid is quite arbitrary and we are not forced to use too small 
time intervals. We find very good agreement between the MC 
method and the analytical predictions. A key point of this algo- 
rithm is that it correctly reproduces the mean statistical quanti- 
ties of EPS and PS predictions for any set of Mi,zi, Mq, Zo- This 
means that there is no need to divide the time interval between zo 
and zi in smaller time-steps. This is particularly convenient be- 
cause it allows us to choose the time grid for the semi-analytical 
implementation in an arbitrary way. 

The algorithm described is based on the Extended Press- 
Schechter formalism, that means on the spherical barrier model, 
but it can be extended for a more general moving barrier model, 
useful for the descrip tion of the ellipsoidal collapse model 
("Bond & Myers', '1996') that bears closer consistency with N- 
body simulation ( Sheth & Tormen, 2002; Gio coli et al.L l2007t 
iZentnerl 120071: IZhang. Ma. & Fakhourii 120081) . 



3. The baryonic model 

While the theory of the hierarchical development of the dark 
matter structures is relatively simple and now well established, 
the complex physics of baryons presents still highly uncertain 
aspects, and requires the treatment of several processes in a sim- 
plified form. The latter is obtained with the introduction of a 
minimum number of free parameters in the model, that can only 
be set through a comparison of the model predictions with ob- 
servational data. These parameters flag our lack of knowledge 
about the details of the physics of galaxy formation. Obviously, 
although they are set by fitting model results with observations, 
only limited ranges of values are allowed: values outside these 
would indicate fundamental problems in the process treatment or 
lack of important physical ingredients. This Section is dedicated 
to a description of our new SAM model, GECO. 
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At variance with the majority of the models of galaxy forma- 
tion, we tested our results directly on the observables involving 
stellar mass rather than luminosities, to avoid introducing further 
uncertainties in the model due to the spectro-photometric syn- 
thesis, dust extinction and the Initial Mass Function. The elab- 
oration of the photometric synthesis based on the model galaxy 
mass function will be deferred to a future work. 

3.1. Linking Galaxies to Dark Matter Halos 

As described in the previous section, the generation of the 
merger tree proceeds backwards in time. Starting from the initial 
redshift (typically z - 0), we split each halo into progenitors at 
higher redshifts. On the contrary, galaxy formation is modelled 
forward in time, starting from the "leaves" of the tree, i.e. halos 
at high redshift whose progenitors fall below the mass limit res- 
olution (Mres = IO'^Mq in our case) and their hierarchy is no 
more followed in time. Starting from each of these halos we put 
baryons inside them, according to the baryonic fra ction observed 
by WMAP5: fb = OJO,,, = 0.17 ('Dunkle5Leta]l|2009). 

We use the grid in parent halo masses and time-steps de- 
scribed in Sect. 2.2. To solve the differential equations for baryon 
evolution that we are going to introduce in the following sec- 
tions, we divide each step of the merger hierarchy into 100 
smaller substeps which allow us to solve the system of equations 
through the use of finite increments of mass. 

In order to derive statistical quantities such as the mass func- 
tion or the SFR density, we simply make a weighted sum over the 
model galaxies. Since we are dealing with an EPS merger tree, 
we use for consistency as a weight scheme the number density 
of halos given by the PS mass function. In some cases, we have 
also checked the results using a Sheth & Tormen number density 
dSheth & Tormeni 12002 ). which fits the results of N-body simu- 
lations with higher accuracy, but we have found that it makes 
little difference. 

We will describe in the following the baryonic processes at 
work in driving the evolution of galaxies, starting from the cool- 
ing of the gas. 

3.2. Dynamical Time-scales and Gas Cooling 

The cooling of the gas is a fundamental process for galaxy for- 
mation, as it sets the rate at which gas becomes available for 
star formation. The first to invoke the need of some dissipative 
process occurring inside dark matter halos were White & Rees 
(1197 8). in order to explain the difference between the luminosity 
function of galaxies, with its characteristic mass and size, and the 
halo mass function, that on galactic scales is a typical power-law 
without any preferred scales. White & Rees argued that galaxy 
formation is a two stage process, with dark matter halos forming 
through the dissipationless hierarchical clustering, and with gas 
cooling occurring inside them. 

The gas is assumed to be heated by shocks during the violent 
relaxation following the collapse, up to the virial temperature of 
the halo, given by; 



2k" 



(6) 



where yu is the mean molecular mass of the gas, mn is the 
mass of a hydrogen atom, k is the Boltzmann's constant 
and Vc is the circular velocity of the halo, that can be re- 
lated to the mass of the halo through the expression: Vc - 
y/GMhaio/Rvir- The virial radius, is the limiting radius 



of a virial halo, within which the mean density is 200pe, 
where pc is the universal critical density, and is given by 
R,,r ^ 1.63 X 10-^(Mw„/M o/2-^'/3(Q(z)/Qj'/3(l +z)-'h-'kpc 
(In avarro. Frenk. & WhiteLll997h . From Vc and the virial radius 
Rvir, an halo dynamical time can be defined: 



thalo - Rvir/Vc ■ 



(7) 



The rate at which the gas cools depends upon this tempera- 
ture, which determines its ionization state, and on the metallicity 
of the gas, that is the chemical comp osition. We take the cool - 
ing function A(T,Z) tabulated bv Sutherland & Dopital (Il993h . 
Note that the cooling function also depends on the metallicity 
of the hot gas. High metallicity values increase the cooling rate, 
because of line cooling from heavy elements, mainly for low- 
mass halos. We assume that the hot gas has a constant metallic- 
ity equal to 0.3Zp , close to the value for the hot gas in cluster 
(ISomerville et al.L 12008). We have also checked that by adopting 
a solar metallicity instead of the subsolar one, the results do not 
change significantly. 

The cooling time is defined as the thermal energy density 
divided by the cooling rate per unit volume: 



^cool 



3 ^ninkT 
2A(r,Z)pg, 



(8) 



where pgas is the gas density profile, which is assumed to be 
isothermal: 



Pgas(r) = 



Mhc 



47TR,i,.r 



2 ' 



(9) 



where M/,o, is the total gas mass in the hot component, which 
is assumed to extend to the virial radius. A cooling radius, rcooi, 
can now be defined as the point where the local cooling time is 
equal to the age of the universe at that epoch. 

We can now compute the rate at which gas accretes to 
the centre of the halo, becomin g available for sta r formation. 
Following the treatment of Whit e & FrenkI (Il99lh . we distin- 
guish two regime of cooling, depending on the value of rcooi- 
The first case corresponds to rcoi > Rvir, the cooling radius ly- 
ing outside the virialized region of the halo. The cooling is so 
rapid that the infalling gas never reaches the hydrostatic equilib- 
rium and the supply of cold gas for star formation is limited by 
the infall rate rather than the cooling rate. We assume that all the 
hot gas in the halo will settle to the centre in a time-scale given 
by the halo dynamical time: 



Mc. 



ool 



R\ir 



(10) 



This regime of cooling was called bv lWhite & FrenkI (Il99lh the 
rapid cooling regime, also known as cold mode. 

The second way of infalling is the static hot halo regime, 
or hot mode, that occurs when the cooling radius lies inside the 
virial radius. In this case, the gas inside rcoi will be pressure- 
supported and will contract quasi-statically toward the centre. 
Cooling will cause a flow of gas toward the centre in a way ex- 
actly analogous to the cooling flow occurring in galaxy clusters. 
A simple expression for the infall rate is given by: 



M, 



cool 



4np 



gas 

cool 



cool 

dt 



(11) 



Note that for small halos (that we find preferentially at early 
times), the cooling radius is greater than the virial radius, and we 
are in the cold mode regime, where equation[TO] applies. At late 



E. Ricciardelli and A. Franceschini: GECO: Galaxy Evolution COde 



5 



times, for large values of the virial velocity, the cooling radius 
falls below the virial radius and we are in the hot mode regime. 
As a consequence, at high redshift we have a more e fficient cool- 
ing. H ydrodynamical simulations, such as that of IKeres et al.l 
(|2005|), justify the use of two different modes of cooling. Indeed, 
they found that a cold mode dominates for low-mass objects, that 
are found at high redshift or in low-density environments in the 
nearby universe, while the hot mode contributes significantly for 
high-mass systems, therefore being important only at low red- 
shift. Anyway it is worth to note that in our model the cooling at 
high redshift might be enhanced due to the assumption of a con- 
stant metallicity. In case of a metal-dependent cooling rate we 
would obtain a lower cooling in the high redshift objects with 
primordial metallicity. We will test such effect in a future ver- 
sion of the code, where the chemical evolution model will be 
implemented. 

3.3. Galaxy Sizes and Angular Momentum 

Since several time-scales involved in the analytic treatment of 
galaxy formation (such as for example the star formation rate) 
that we are going to introduce in the following sections depend 
on the dynamical time of the galactic disc (fjv« - rdisc Iodise), it 
is important to have an accurate description of their sizes. The 
disc size will depend on the virial radius of the halo where the 
galaxy was born, and on its angular momentum. Indeed, if the 
halo is asymmetric and surrounded by a clumpy distribution of 
matter, then it can acquire an angular momentum from a net tidal 
torque. To quantify the angular momentum of the system, one of- 
ten refers to the spin parameter, which is a dimensionless quan- 
tity, defined as: 



A = 



J\E\ 



1/2 



GM\ 



5/2 
halo 



(12) 



where J, E and Mhah are the total angular momentum, gravita- 
tional binding energy and mass of the halo. The distrib ution of 
the sp i n parameter as found in N-body simulations ( Warren et alj 
Il992t ICole&Lacevlll996h can be approximated with a log- 
normal distribution with < A >- 0.05 and cr; = 0.5. Hereafter 
we will assume the mean value for A. Mo, Mao, & White ( 1998) 
related the disc radius to the spin parameter and to the virial ra- 
dius assuming that the angular momentum of the disc is a fixed 
fraction of that of the halo: 



- M-y 

^/2\mdJ 



\AR„ 



(13) 



where nid is the ratio between the disc mass and the halo mass, 
and jd is the ratio between the angular momentum of the disc 
and of the halo. Higher values for the spin parameter result in 
larger discs, because they contract less before reaching the cen- 
trifugal equilibrium. The above relation holds in the case of halos 
approximated with isothermal profile and neglecting the gravita- 
tional effects of the discs themselves, but the authors also give 
the correct expression in the c ase of NFW profiles an d self- 
gravitating discs. As shown by Mo. Mao. & Whit3(ll998h . it is a 
reasonable assumption to take jd - nid, which gives us: 



A 

^ disc ~ "^t^vir •> 

V2 



(14) 



which leads to galactic discs that are over an order of magnitude 
less extended than their halo hosts. Under the assumption that 
we are considering negligible the self-gravity of the disc, the 
disc circular velocity is equal to the halo virial velocity. 



The assumption that the specific angular momentum of the 
disc is the same of that of the halo, as implied by jd - 
wid, leads to disc s i zes th at well match the observed spirals 
dM o. Mao. & Whitel Il998h . hence it became the standard as- 
sumption of disc modelling. Nevertheless, hydrodynamical sim- 
ulations find that the gas looses most of its angular momentum 
during accretion, producing discs too small compared to that ob- 
served (Navarro & White, 1994). It is still unclear if this effect 
is simply due to the limited mass resolution employed, or to the 
absence of some form of feedback. 



3.4. Star Formation 

Obviously, the physical processes governing the rate with which 
the (primordial) gas is transformed into stars play a critical role 
in models of galaxy formation. Unfortunately, they are still sub- 
ject to major uncertainties, given the extremely complex inter- 
play of phenomena on the largest (magnetic fields, gas angular 
momentum, etc.), as well as smallest (radiative cooling, thermal 
conduction, dust extinction, etc.) scales. Nevertheless, observa- 
tions indicate the existence of some useful laws to describe SF 
on global scale, larger than star-cluster scales. Anal jysing a large 
sample of spirals and starburst galaxies, Kennicutt' (1998) found 
two empirical correlations for the star formation per unit area, 
I.SFR- (i) ^SFR ^d (ii) I.s FR ^gas/tdyn- Thcsc empir- 

ical estimates justify the use of parametric forms for the star- 
formation prescription in semi-analytical models, where the star 
formation time-scale is given by the dynamical time. 

In our model we consider two modes of star formation. In 
the first one, the quiescent mode, the star formation rate is sim- 
ply assumed to be proportional to the amount of cold gas avail- 
able and inversely proportional to the time-scale of star forma- 
tion, that is assumed here to be the dynamical time of the disc: 
td\n - r disc I V disc- The efficiency of the star formation process is 
quantified by the free parameter aquie defined by the relation: 



M, 



cool 



^dvn 



(15) 



Since this gas is assumed to share the same high value of the 
angular momentum as previously discussed, stars formed in this 
way are added to the disc component of the galaxy. 

The second mode of star formation considered in GECO is 
the burst mode. Starbursts are allowed to occur after a merger 
between galaxies, that can be due to the cannibalism of a satel- 
lite galaxy by the central one or to a collision between satel- 
lites. Indeed, a strong correlation between galaxy interactions 
and starburst activity is observed both in the present-day uni- 
verse and at high redshift. This is for example the case for 
ULIRGs, the Ultra-Luminous Infrared Galaxies discovered lo- 
cally by IRAS and then found in large numbers at high-z by ISO 
and Spitzer, exhibiting the most violent ep isodes of star forma- 
tion (.Kennicutt, ,1998; .Sanders & Mirabel! J 996,) . The assump- 
tion of a starburst triggered by a galaxy m erger is furthermore 
justified by several numerical simulations (iMihos & HernquistL 
1994, 1996), that show how, during the merger phase, gas 
sinks to the central region of the galaxy and the increased den- 
sity leads to an enhanced level of star fo rmation (accordin g to 
the Kennicutt law). Recent simulations dCox et al.L l2008h in- 
dicate that the starburst event can be triggered in the major, 
as well in the minor, mergers, but with an efficiency which 
strongl y depends on the galaxy mass rat i o (as implemente d 
also in ISomerville. Primack. & Fabeiil2001l:rMenci et alJl20"04l) . 
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Therefore, in our model, the star formation during the starburst 
is given by: 



burst ~ 



^burst 



(16) 



where thursi is the dynamical time of the largest disc, while f burst 
is the efficiency of the burst, assumed to be linearly proportional 
to the merger mass ratio: 



f burst — P burst ? 

m2 



(17) 



where nii and m2 are, respectively, the smallest and the biggest 
galaxy, and /3burst is a free parameter. Stars formed in such way 
are added to the bulge component of the galaxy. 

In addition to the increased rate of star formation, major 
mergers between spirals move stars from circular orbits to ran- 
dom motions, hence producing a remnant that generally resem- 
bles an elliptical galaxy. For these reasons, we can safely assume 
that after that kind of mergers, the discs previously present in 
the colliding systems are transformed in a spheroid. We define a 
merger to be a major one when > 1/3. Spheroid formation 
through disc instability is not considered in the present work. 

3.5. Feedback Processes 

Feedback processes are mechanisms which regulate the effi- 
ciency of star formation. The need for some form of feedback in 
low-mass hal os, which suppresses c old gas accretion, was first 
recognised bv lWhite&ReesI (119781) . in order to flatten in some 
way the faint-end of the luminosity function, which turned out 
to be too steep compared to the observational data available at 
the time. SN feedback and photoionization are commonly con- 
sidered to reduce the faint-end slope. 

Although initially the motivation for invoking feedback was 
to suppress cooling, and hence star formation, in low-mass ha- 
los, in recent years the focus was shifted to attempt reproduc- 
ing the bright-end of the luminosity function, which in the most 
classical models resulted in an excess of bright objects at low 
redshift compared with the observations (again a manifestation 
of the basic problem to shape an essentially power-law form for 
the DM mass function into a Schechter exponential function for 
the baryons). Although the behaviour of the cooling function im- 
plies that the cooling time increases in high-mass halos, this is 
not enough to explain the sharp observed cut-off in the galaxy 
luminosity or stellar mass functions, and some mechanisms pre- 
venting the cooling rate in high-mass halos is also needed (and 
is commonly identified in the AGN effect, see below). 

3.5.1. SN feedback 

The first form of feedback considered in our model is that from 
SN explosions and high-mass stars outflows, which eject gas 
and energy into the surrounding Interstellar Medium (ISM). This 
is supported by several eviden ces of Supernova driven winds 
in dwarf star-forming galaxies (Martini 119991: [Strickland et al.L 
[2000; Heckman, 2001), which also suggests that the gas reheat- 
ing rate is proportional to the SFR. 

The parametrisation of the SN feedback is based on sim- 
ple energy arguments. The rate of reheating will be propor- 
tional to the SFR, to the number of Supernovae per solar mass 
of stars, given by rjsN = 4 x for a Salpeter Initial 

Mass Function (IMF), to the energy released by each SN, about 
EsN = IQ^^ergs, and to the efficiency with which the SN energy 



is transferred into the ISM, esN- This last parameter is highly 
uncertain, and is usually treated as a free parameter The rate of 
reheating will be more efficient for galaxies living in low poten- 
tial wells. In order to mark the depths of the potential well of the 
host halo, we use the halo circular velocity, so that the rate of 
reheating will be simply proportional to V^^. This implies that 
galaxies living in low-mass halos are affected by SN feedback 
effects, and their star formation will be self-regulated, while in 
high-mass halos this kind of feedback is ineffective. The rate of 
reheating will be finally given by: 



JyJheal - eSN — y2 



(18) 
[l99l 



(Kauffmann, White, & Guid erdonil 
Somerville & Primach i 1999r 

The reheated gas is removed from the cold gas reservoir, and 
can not be used to form new stars. We assume an ejection model 
of feedback, that means that the reheated gas is ejected from 
the halo and is unable to cool until it is reincorporated into a 
more massive halo at the following step of the merging hierar- 
chy. Although our re-incorporation scheme assurnes a fast fall- 
back (analogous to the fast re-incorporation of iDe Lucia et al.l 
|2004), we found that it produces better results compared to the 
retention model, in which the gas can subsequently cool. In par- 
ticular we obtain a better match of the faint-end of the stellar 
mass function. 

3.5.2. Reionization 

It is now known that the hydrogen in the intergalactic medium 
(IGM), that became neutral at z ~ 1000, must have been ionized 
at later epochs, although the redshift at which this reionization 
occurred is still quite uncertain, ranging from z = 6, as imposed 
by the lack of a Gunn-Peterson trough in quasar spectra at that 
redshift ( Fan et al.L 12000 ). up to z = 1 1, as imposed by the con- 
straint on the optical depth to the last scattering surface measured 
from WMAP5 data. 

If a large population of galaxies and quasars exist at very 
high redshift, as predicted by galax y formation mod els and con- 
firmed up to z ~ 6 by observations dFan et alil2000l) . then reion- 
ization could have occurred through photoionization, since both 
young galaxies and quasars emit UV photons, able to ionize 
the IGM. This photoionizing background may also act to inhibit 
galaxy formation in two ways. It can heat the IGM, increasing 
the thermal pressure of the gas and preventing it to collapse into 
the dark matter halos. The second way is through a reduction 
of the cooling rate of gas inside halos, mainly reducing neutral 
atoms which can be collisionally excited. This results in a strong 
suppression of galaxy formation in l ow-mass halos. 

According to previous studi es (iThoul & Weinbergl 119961: 
iGnedinL l200a [S omervi iia i2ool . after reionization the above 
mechanisms prevent gas accretion in halos with Vc < 50 km s \ 
with the result that galaxy formation is strongly suppressed in 
such systems. Hence, in our model, we mimic the effect of pho- 
toionization by suppressing gas cooling in halos having a cir- 
cular velocity lower than that at redshifts greater than the red- 
shift of reionization, assumed to be Zre - 6. Although the exact 
value of the redshift at which reionization starts is still uncer- 
tain and other models treat it as a free parameter, we realized 
that by varying its value our results does not change apprecia- 
bly and we choose to fix it. This is in agreement with previ- 
ous works that found that the properties of the galaxy popula- 
tion are almost insensitive to the assumed redshift of r e ioniza - 
tion dKravtsov. Gnedin. & KlvpinL 12004 iMaccio et all l2009l) . 
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Indeed, at very high redshift the leaves of the tree have very 
small mass and the only cooling mechanism is the H2 cooling, 
which is very inefficient. 

The simple approa ch used in ou r mode l for the reionization 
were demonstrated by Benson et alj (l2003b to have very similar 
effects to the self-consistent treatment of the photoionization de- 
Benson et al. (2002). In Appendix B we 
present the effect of varying the assumption on the limit of the 
circular velocity. 

3.5.3. AGN feedback 

There is growing evidence for a tight relationship between 
galaxy evolution and the growth of super-massive black holes 
(SMBHs) powering nuclear activity. Supporting evidences come 
from the tight co rrelations betwee n the BH mass and the mass 
of the bulge (Har ing & RixL |2004 and that between the black 
hole mass and the velocity dispers ion of stars in the bulge 
dFerrarese & Merrittl IJOOO; Gebhard t"etall l2000h . The mutual 
feedback between galaxies and quasars may be the reason for 
such strong correlations. 

We implement the accretio n onto black holes in our model 
foll owing the pr escriptions of iKauffmann & Haehnelti (l2000l) 
and lCroton et al.l (2006). We allow two different modes of feed- 
back, one following the starburst triggered by galaxy mergers 
and strong interactions, the quasar mode, and a second one, the 
radio mode, taking place at low rate during the whole life of 
the galaxy. A third natural way of growth for SMBH's is repre- 
sented by the coalescence of the BHs residing in the centre of 
two merging galaxies. In this case, after the merger of the host 
galaxies, the new BH mass is simply the sum of the two progen- 
itors. We thus simplify our treatment by ignoring the occurrence 
of BH binary systems and gravitational wave losses of angular 
momentum. 

We use as a starting point a seed mass for the BH in the 
galaxies on the bottom of hierarchy, i.e. living in the "leaves" 
halos. These seeds are believed to form at very high redshift, 
but it is quite unclear if they formed from the direct collapse of 
pre-galactic gaseous discs or if they are the remnants of ma ssive 
Population III stars dVolonteri. Lodato. & Nataraianll2008l) . We 
take a seed mass equal to IOOOMq, and checked the consequence 
of this choice by using different mass values: the results do not 
depend on this choice, at least for small values of the seed. Since 
we start placing the seeds at the bottom of the hierarchy, the 
redshift at wich they are planted depends upon the time at wich 
the leave haloes fall below the mass resolution of the tree and it 
can be higher than the redshift of the reionization. Anyway, they 
can not grow before the reionization starts since gas accretion, 
that is the main mechanism of feeding BHs, is inhibited. 

The first mode of accretion is motivated by the observations 
of bright AGN and quasars radiating with high efficiency, close 
to their Eddington limit. It is widely believed that this mecha- 
nism of accretion is triggered by galaxy mergers, since they can 
drive rapid inflow of gas toward the centre of the galaxy, feeding 
both the starburst and the BH. 

As in 'Kauffmann & Haehnelti (I2000'). we assume that during 
a merger of galaxies a certain fraction of cold gas is accreted on 
to the centre of the black hole: 



QSO 



■old 



1 -H(280kms-VVc)2 



-At" 



(19) 



where fsH is an efficiency parameter to be chosen in order to 
match the relation between black-hole mass and velocity disper- 
sion of the bulges and At is the time interval of the integration 



of differential equations. In this way QSO accretion is closely 
linked to the starburst activity of the galaxy, which is triggered 
by galaxy mergers together with BH accretion, naturally pro- 
ducing BH masses proportional to the mass of the spheroidal 
component that is formed only during the starburst. 

The second way of accretion, the radio mode, is a continuous 
and quiescent accretion during the whole life of the galaxy, dur- 
ing which the black-hole accretes gas directly from the hot halo. 
This is supported by observations of radio galaxies (mainly in 
high-density environments, like groups or clusters of galaxies), 
which are accreting at a sub-Eddington rate, hence in a rather 
inefficient fashion. In spite of this inefficiency, the energy ex- 
tracted from the BH accretion in the form of collimated jets is 
believed to contrast the cooling flow in the centre of clusters. 
Assuming that in this case the BH is fuelled with a Bond i accre- 
tion rate, the resulting BH gro wth rate is given by (.Croton et all 
lloOetlSomerville et alil2008h : 



RADIO 
BH 



fhotl Vc \^ 

.l\200kms-i/ ' 



108Mo0.1\200kms 



(20) 



where is the ratio between the hot gas mass and the halo 
mass and kACN is a free parameter, with the dimensions of an 
accretion rate (Mq yr"'). 

The most important channel building up the mass of the 
black hole is the quasar mode, in which the accretion rate can 
be much greater than the Eddington luminosity. The second way 
of accretion, the radio mode, provides negligible contribution to 
the present-day black-hole masses, being several order of magni- 
tude below the Eddington rate. Nevertheless this is an important 
source of feedback in high mass halos, as it suppresses cooling 
flows as well as the quasar mode mechanism, and it occurs for 
whole life of the galaxy, while the quasar mode acts only in a 
very short period of the galaxy life. The feedback efficiencies of 
such processes, very difficult to model a-priori, are set by trying 
to match various observables. The efficiency of the quasar mode 
accretion is set by requiring the model to match the MgH-Mhuige 
relation, while the radio mode accretion efficiency is tuned to 
optimize the shape of the stellar mass function, and in particular 
the observed knee at the characteristic stellar masses. 

The injection of energy into the ISM due to the presence 
of an AGN is modelled in the following way. The mechanical 
heating generated by the the black hole accretion is given by 



Lbh = rfMsHC 



(21) 



where c is the speed of light and 77 is the global efficiency of 
energy production close to the event horizon in nuclear black- 
holes, and is set to the reference value of 0.1. The gas cooling 
rate in the halo is then corrected by this energy injection accord- 
ing to: 



M' 



old 
cool 



(22) 



To avoid an unphysical negative value for the cold gas accre- 
tion rate, we enforce mass conservation by assuming an amount 
of newly formed cold gas mass in the time-step At the max- 
imum between zero and M'.'™^At, where M"™^ is given by eq. 
(I22I 1. Hence, the cooling rate can not only be reduced, but even 
stopped in the case of a strong BH accretion rate. 

3.6. Galaxy Dynamics 

During the merger of two dark matter halos, the baryonic cores 
that they contain, being more compact and less subject to tidal 
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effects, may avoid merging with each other, and end up orbiting 
within the new combined halo. An halo formed by many merg- 
ers may contain many distinct galaxies, with the galaxy in the 
most massive halo becoming the central galaxy and all the others 
becoming the satellites. These satellite galaxies gradually loose 
their energy and angular momentum under the action of dynam- 
ical friction until they sink to the centre of the halo and merge 
with the central galaxy. Moreover, they are subject to random 
gravitational collisions that can occur with other satellites. 

3.6.1. Dynamical Friction 

The effect of dynamical friction in slowing down the orbital ve- 
locities of galaxies moving throug h a sea of dark matter par- 
ticles was originally described by ' Chandrasekhad (1194 3*) and 
previous SA M made use of this approach. A lthough early N- 
body studies (iNavarro. Frenk. & Whitel Il995h gave support to 
the merger time-scales predicted by the classical appro ach, re- 
cent numerical and hydrodynamical s imulations (Jian g et al.L 
I200 8': Boylan-Kolchin, Ma, & Ouataert, 2008) have shown that 
the Chandrasekhar's approach underestimates the merger time- 
scale for minor mergers and overestimates it for major merger 
events. As an estimate of the dyna mical friction time scale we 
have assumed the fitting formula of iJiang et al.l (l2008l) . derived 
from hydro/N-body simulation: 



tdf = 



0.94e" *' + 0.6 Mh 



halo 



1 R., 



2C 



m ln(A) Vc 



(23) 



where C is a constant equal to 0.43, In A is the Coulomb log- 
arithm, that can be approximated to In A = ln(l -i- M/,„/o/m), 
e is the circularity parameter of the satellite's orbit, given by 
e = J/Jc, with J the angular momentum of the actual orbit and 
Jc the one of the circular orbit with the same energy. We adopt 
the a verage value found in numerical simulations: < e >= 0.51 
d^men, 1997; Jiang et al., 2008). 

It is worth noting that the time-scale increases as the halo 
mass increases, so that at high redshift, when the halos are less 
massive, we find an high rate of merging, whereas it decreases at 
low redshift, where we have an accumulation of sub-halos inside 
the parent halo. 

An important consideration concerns the mass m used to 
compute the dynamical friction time-scale. As soon as a galaxy 
becomes a satellite in a larger halo, its initial mass is deter- 
mined by its total (dark matter and baryonic) mass. However, 
as it shrinks towards the centre, a large fraction of the dark halo 
is stripped away, while the mass of the gaseous core is too dense 
to be stripped. Thus, the effective mass of the satellite decreases 
with time and the dynamical friction time-scale should increase. 
We realized that, by neglecting this effect, the time-scales for 
merging were too small, leading to an overmerging, particularly 
at low redshifts. We have taken into account the effect of tidal 
stripping allowing the mass of satellite varying with time accord- 
ing to the time-scales predicted by Taffpni^lal. (2003). Hence 
we assume that the satellite looses their dark halo 1 .5-2 Gyr since 
entering the virial radius (for simplicity here we are neglecting 
the tidal stripping on the baryonic core, hence we do not consider 
tidally destroyed satellites). 

3.6.2. Satellite Collisions 

In addition to the mergers driven by dynamical friction, we in- 
clude random collisions between satel htes. We consider for the 
satellite encounters the cross-section of lMamo ^ (Il992h . which is 



a simplified version of t hat derived from N-body simulations by 
IRoos & Normanl(ll979h . The resulting merger rate is given by: 



k = 2 ^f^ala,Rl,.VgalK{yJvgal) , 



(24) 



where Vgai is the internal dispersion of the satellite galaxy, Vc and 
Rvir are the circular velocity and virial radius of the host halo, 
respectively, ap - 4 and a,, = 5.4 are dimensionless parameters, 
while K isa dimensionless merger rate given by: 



K 



1 



- + — 



2 exp(- 



-x2) 



3 -\/n erf(x) 



OvVgal 

2Vc ' 



(25) 



(26) 



For realistic values of the ratio Vc/vgai the merger rate decreases 
with this ratio and becomes very low for galaxy clusters. Given 
the merger rate we can compute the collision time as: 



tcoll 



1 

kn 



(27) 



where n is the number densities of galaxies in the halo. We ob- 
tain shorter merger time-scale at high redshift, when the num- 
ber density of galaxies is higher. The probability for having a 
merger for a satellite with another one, randomly chosen, will 
be P = At /tcoll, where At is the lifetime of the halo. We consider 
only binary mergers, and we assume that such event occurs when 
f > 1 . During a collision between two satellites we allow a burst 
of star formation in the same way as the one occurring during a 
merger between a satellite and the central galaxy. The effect of 
including such process, in addition to the mergers involving the 
central galaxy due to dynamical friction, is to decrease the frac- 
tion of low-mass galaxies in favour of intermediate-high mass 
objects, hence to slightly modify the shape of the galaxy mass 
function. 

3.6.3. Remnant 

After the merger, the total mass of the remnant (dark matter and 
baryons) is the sum of the two merging galaxies, while the new 
radius and circular velocities are computed from applying the 
conservation of energy and the virial theorem. According to the 
virial theorem, the total internal energy is given by E,,,, - -T. 
By applying the conservation of energy: 



-T\ +T2- Eorh , 



(28) 



where T denotes the kinetic energy and Eo,-b is the mutual orbital 
energy: 



Eorb — forb 



Gmim2 
r\ + r2 



(29) 



(fx and r2 are the radius of the two progenitors and fort is a pa- 
rameter which weakly depends on the density profile, we assume 
forb - 2). These considerations yield to: 



Rncw — 



(nil + OT2) 



myri + m\lr2 + forbmimliri + ^2) 



(30) 



In the case of major mergers involving two equal galaxies (nii - 
m2 and ri - r2) the size of the remnant is a factor of 4/3 greater 
than the initial size. 



Table 1. Free parameters 
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Parameter Meaning Best-fit value Range 

Q-^,„,> star formation 0.01 0.01-0.1 

Pburst burst efficiency 0.68 0.5-1 

SN feedback efficiency 0.5 0.01-1 

fBH AGN feedback efficiency in QSO mode 0.006 0.001-0.1 

kAGN AGN feedback efficiency in radio mode fx 10"* 10"^ - 10"* 




Fig. 2. Left panel: Stellar mass function at z-Q. Black solid line represents the SMF resulting from GECO, purple shaded region 
stays for the COl observed mass function with its error-bar, and the green shaded one from B03. Also shown is the contribution 
to the total stellar mass content from bulges (red long dashed line) and discs (blue short dashed line). Right panel: Stellar mass 
function at z-0 of early-type (red long dashed line) and late-type galaxies (blue short dashed line) compared with the observed one 
by B03: red squares for early-type and blue triangles for late-type. 



3. 7. Galaxy Morphology 



In our model the morphology of a galaxy is determined by the 
relative importance of the bulge over the disc component. The 
only way of formation for discs is through the quiescent mode of 
star formation, while bulges can grow in two ways: through star 
formation occurring in a starburst event, hence triggered by any 
event of merger, and through the disruption of discs following a 
major merger In this case all the stars belonging to the discs un- 
dergoing a merger are added to the bulge of the remnant galaxy. 
Such galaxy may eventually form a new disc if some fraction of 
cold gas is still present. 



We can assign a crude morphological type to each galaxy 
by using the ratio between the bulge mass and the total 
stellar mass: r - MtuhelMstar- U sing the prescription of 
iBertone. De Lucia. & Thomas! (l2007h . we classify as Ellipticals 
galaxies with more than 70% of their stars in bulge, as Spirals 
galaxies having < r < 0.7, and as Irregulars galaxies without 
any bulge. 



4. Setting the free parameters 

In the previous section we have introduced several free param- 
eters in the analytical prescriptions which model galaxy forma- 
tion. In summary they are: 

1. aqiiie'. efficiency of star formation in quiescent mode 

2. Pbursi'- efficiency of star formation in starbursts 

3. esN'- efficiency of SN feedback 

4. fsH'- efficiency of AGN feedback in the quasar mode 

5. kAGN'- efficiency of AGN feedback in the radio mode 

We will first infer constraints on their values by comparing the 
model results, obtained for a large grid of parameter values, with 
observations of the local universe. In particular we focus on the 
build-up of the stellar mass, and we require our model to fit the 
local stellar mass function, as well as the relationship between 
the black-hole mass and the bulge mass. The stellar mass func- 
tion is obviously influenced by all the parameters together The 
star formation efficiency in the quiescent mode Oquie affects the 
overall shape and normalisation of the function; esN affects the 
faint-end (or low-mass end); fibursi determines the contribution 
of the bulge on the total stellar content, and we require that at 
the bright-end bulges dominate over discs. The AGN efficien- 
cies, fsH and kAGN shape the bright-end of the mass function. 
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Fig. 3. Black-hole versus bulge mass relation at z-Q. Blue tri- 
angles represent data for modelled galaxies, while black squares 
with error-bars are data from .Haring & Rix. (2004.) . The magenta 
solid Une is their best fit relation. 



Moreover, the AGN efficiency in the quasar mode, i.e. the one 
triggered by the starburst episode, determines the relationship 
between the mass of the black-hole and the bulge mass, affecting 
mainly the normalisation of the relation, while the slope turns 
out to be very close to the one observed almost independently 
on the value of this efficiency, since in the model the growth of 
BHs and bulges are closely linked. In Table[T]we report the value 
of these parameters in our best-fit model and the allowed range. 
Unless otherwise stated, in the following we always will refer to 
this best-fit model. 



5. Matching data from the local universe 

A key observable for comparing predictions of any models of 
galaxy formation is the number density of galaxies as a func- 
tion of their stellar mass in the local universe, the stellar mass 
function (SMF) at z = 0. We will use it to constrain our model 
free parameters. We refer to two measurements o f the local 
SMF. The first is that reported by ICole et al.1 (1200 ih . hereafter 
COl, obtained by combining data from the Two Micron All Sky 
Survey (2MASS) and the 2dF G alaxy Redshift Survey. The sec- 
ond estimate, by Bell et al. (2003), exploited a large sample from 
2MASS and the Sloan Digital Sky Survey (SDSS). We choose 
to show their ^-selected sample (taken from their table 5). Both 
estimates use near-infrared data to compute stellar masses and 
hence are highly reliable. All data are transformed to a Salpeter 
IMF 

In the left panel of Figure |2] we show the stellar mass func- 
tion at z = 0, resulting from the model and compared with 
the COl and B03 observations. The model is able to reproduce 
with good accuracy the shape and normalisation of the SMF, 
especially in the bright-end {M^tar > IO^'Mq/i"'). Compared 
to the COl mass function, the model shows a slight excess of 




Fig. 4. Local Black-hole mass functio n of the model ( 'black solid 
line) compared with the estimates of IShankar et alJ (l2004l) (red 
shaded region). 



galaxies at the faint-end. However, the BOB estimate presents a 
steeper faint-end, resulting in closer agreement with our predic- 
tion. Our model also shows an underestimation of galaxies at 
the knee of the mass function (lO"'-^ - IO'^Mq/Zi^), an eff^ect 
already present, for instance, with the Mille nnium galaxy cat- 
alogue ( Bertone. De Lucia. & Thomasl l2007l) . Such dearth may 
be partly explained with the combined effect of SN and AGN 
feedback that are both effective in this mass range. Moreover, 
since the deficit is due to a lack of intermediate-mass early-type 
galaxies, (see right panel of Fig. |2]i it is likely connected to the 
excess of low-mass spheroids, as satellite galaxies are not able to 
become massive enough in the model for the reasons described 
later in this section. 

It is worthwhile to stress the fundamentally different shapes 
of the halo mass and the luminous (stellar) mass functions pre- 
sented here. The general shape of the dark matter halo mass 
function (see Figure [U, as predicted by any theoretical models 
based on hierarchical dark-matter clustering, is a pure power- 
law with a steep faint-end and no knee at the high masses. A 
turn-off mass does exist, but it occurs at much higher masses 
(M ~ 10'^ Mq) than the characteristic one observed in the stel- 
lar mass function. In order to obtain a Schechter shape for the 
mass function, many effects contribute. At the faint-end side, 
the feedback from SN and from the photoionization background 
suppresses star formation in small halos, hence reducing the 
faint-end slope. At the bright-end, the difference between the 
exponential cut-off observed in the stellar mass function and 
the power-law shape expected for the halo mass function is ex- 
plained with the inter-play of various independent effects. On 
one hand, there is a marked dependence of the gas cooling on 
the halo mass: the cooling time decreases with the halo mass. 
In addition, the two modes of gas cooling discussed in Sect. 4.3 
have an halo mass dependence: in the cold mode regime the gas 
accretes towards the centre in a time-scale given by the free-fall 
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time, while in the hot mode regime the accretion rate is governed 
by the cooling time, and is much less efficient than in the for- 
mer case. Since in high-mass halos the accretion mainly occurs 
through this last mechanism, this may explain why in this halos 
gas is not able to cool efficiently. Moreover, the AGN feedback 
acting on this scales, further reduces the cooling, leading to an 
approximately correct number of massive galaxies. 

In Figure |2] (left panel) the separate contributions to the 
total stellar content from discs and bulges are also shown, as 
blue short-dashed and red long-dashed lines, respectively. As 
expected, the bright-end is dominated by bulges, meaning that 
objects with very high-mass {Mstar > IO^'Mq/z"') do not essen- 
tially contain discs, that have been presumably destroyed dur- 
ing a major merger. On the contrary, discs overcome the num- 
ber densities of bulges at low masses, and start to be domi- 
nant in the mass range Ms,„r ~ IO^'-^Mq/i"'. Note that this 
value is close to the transition mass, ~ 3 x IO^^Mq/;"' at 
z-Q, at which a transition in s everal galaxy properties is ob- 
served dKauffmann et all l2003h and where early-type galaxies 
start to d ominate over the late-type ones in the local universe 
(iBaldrv et al., 2004; Bundy, ElHs, & Conselice, 2005]). 

So far we have considered only the contribution to the total 
stellar mass density from the disc and bulge components. In the 
right-end panel of Fig.|2]we show the mass function of early-type 
and late-type galaxies. As defined in the previous Sect., early- 
type galaxies are defined as galaxies whose bulge makes at least 
70% of the total stellar mass. We include in the definition of late- 
type systems both spiral and irregular galaxies. We compare the 
predicted mass function with the results of B03. With this defini- 
tion of morphologies, the agreement with the observational data 
is quite good, in particular we accurately reproduce the transition 
mass at which the number density of ellipticals/SOs equals that 
of late-type galaxies locally. Anyway, we note quite a large ex- 
cess of faint early-type galaxies with masses below 10"^ Mq/Zi^. 
This is a long-stan ding problem that afflicts all the SAMs in 
the recent literature ([Wei nman n et aL . 2006t Croton et all 120061: 
ICattaneo etall 120081: ISomervill e et all l2008l) and may be re- 
garded as an indication of a wrong treatment of the satellite 
population. A common way among SAMs of treating the satel- 
lites is to assume that the hot gas in a system that merges into 
a more massive one is instantaneously stripped (strangulation) 
and added to the reservoir of the new central galaxy. As a con- 
sequence, the star formation in the new satellite is quenched 
shortly once it has consumed its cold gas. This crude treatment 
of the strangulation is clearly an oversimplification and recent 
hydrodynamical simulations indicate that such a pr ocess occurs 
on a l arger time-scale, of the order of 1-10 Gyr ( Mc Carthv et al.L 
l2008h . allowing the star formation to continue for some time af- 
ter the coalescence, hence lowe ring the fraction of red satellites 
dKang & van den BoschL 2008|). Althou gh even tidal disruption 
(see for instance ISomerville et al.ll2008l) can help in decreasing 
the faint-end, it is unlikely that it ca n change the relativ e num- 
ber of early and late-type satellites ( 'Kimm et al.L |2009). Other 
effects, such as the r am pressure on the cold gas in the galaxy 
dLanzoni et aUl2005l) may also play a role in the formation of the 
satellite population and should be included in the future SAMs. 

In order to set the efflciency of black-hole growth, we require 
our model to match the observed relation between the black-hole 
mass and the mass of the stellar bulge in the local universe. Since 
the major channel for the growth of black-hole is given by the ac- 
cretion during the quasar mode triggered during a starbust event, 
where all the stars formed are added to the bulge component, 
it comes with no surprise that black-hole and bulge are closely 
linked. Nevertheless, in order to obtain an acceptable fit, where 



both the slope and the normalisation of the relation are well re- 
produced, some fine tuning of the /bh and kACN parameters were 
required. 

In Figure [3] we show the comparison between our derived 
relation (blue triangles) with data from Haring & Rix (2004), 
represented by the black squares with error-bars, and their best- 
fit relation, shown by the magenta line. The observed relation 
is well reproduced by our model over the w hole mass range 
(note in any case that the lHaring & Rixl d2004l) fit is derived only 
for bulge masses Mh„igg > IO^^Mq). Anyway, the scatter in the 
model results to be smaller than that observed, as the model ne- 
glects the contribution to the bulge formation due to disk insta- 
bilities, and the only channel to build spheroids is through merg- 
ers, hence directly linked to the BH formation. 

A further test performed on the black-hole population is 
shown in Figure |4] where we show the BH mass function of 
the model, compared with the local estimate bv IShankar et al.l 
(2004). The remarkable agreement between the observed and 
predicted BH number densities, together with the well matched 
BH-bulge relation is encouraging and suggests that our method 
of implementing AGN feedback, though simple, is accurate 
enough. 



6. Comparing GECO with high-redshift galaxy 
statistics 

Having obtained in the previous Sect, a good match of GECO's 
predictions with the global properties of galaxies in the local 
universe, it is instructive now to attempt a first comparison of 
the evolutionary properties of the model with the galaxy statis- 
tics in the distant universe, with no further adjustment of the 
model parameters. Given the wide range of physical processes 
included in GECO, one would expect at least a rough agree- 
ment of the global properties, unless some key physical pro- 
cesses might have been overlooked. 

6. 1. The Evolution of the Stellar Mass Function 

We first test the resulting mass function with observations up 
to z ~ 3.5. A comparison of the GECO's mass functions 
with those observed is reported in Figure [51 Here we match 
our model mass functions with results from Franc eschini et al.l 
(2006), herea fter FR 06, Fontana et al. ( 200_6), hereaf t er F06 , 
Perez-Gonzal ez et al l d2008l) . PG08, iMarchesini et al.l d2009l) . 
M08 and Berta et al.l d2007h . B07. All data are transformed to a 
Salpeter IMF. We consider only mass function estimates which 
take into account IRAC bands in the stellar mass derivation, 
since the lack of Spitzer data in the near-IR leads to an overesti- 
mation of the stellar masses at z > 3 (see F06). Note that error- 
bars in FR06 take into account only Poisson statistics, hence 
resulting to be much smaller than in the other data, where un- 
certainties in photometric redshifts, stellar masses and cosmic 
variance are considered. 

Considering the spread between different estimates of the 
mass function, the comoving number density of massive galax- 
ies in our model is roughly consistent with the observed val- 
ues, although with a tendency to underpredict the data, espe- 
cially at high redshifts. In the highest redshift bin considered 
(z ~ 3.5), the bright-end is marginally consistent with that of 
F06 but systematically lower than M08 and PG07. Note how- 
ever, that including the effects of different assumptions in the 
Spectral Energy Distribution (SED) modelling leads to different 
MF estimates fully consistent with each other (see M08). 
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Fig. 5. Comparison between GECO's mass functions at high redshift, up to z ~ 3.5. The green shaded region indicates the mass func- 
tion o f lpra nceschini et al. 2006, the magent a one is that of Perez-Gonzalez et al. 2008, the blue region is the estimate of jFontana et al] 
l2006l the yellow shading is from Marches ini et al ]^09 and the red squares at z ~ 2 show the results from lBerta et al.ll2b07l GECO 
predictions are represented with a solid black line. The dotted line in all the panels is the GECO mass function at z = 0, reported as 
reference and the red dashed line is the model prediction for the WMAP5 cosmology. 



A remarkable feature in the comparison of model expecta- 
tions with the data concerns the faint end of the mass function, 
where the model tends to keep systematically in excess of the 
observed number of low-mass objects at any redshifts: whereas 
there is a significant evolution in the number density of low 
mass galaxies in the data, the model predicts modest or no evo- 
lution up to redshift of one. The mismatch between data and 
hierarchical models at th e faint-end seems a common feature 
among different models (" Fontana et al.L l2006t iFontanot et al.L 
12009; .Monaco et al...2009.) . 



In the same figure |5] we also include the predictions for 
the SMF adopting the WMAP5 set of parameters indicated by 
the red dashed line. The main diff'erence with the Concordance 
cosmology is a lower value of erg and a tilt in the primordial 
power spectrum (n < 1), implying a delay in the formation of 
structures. Indeed, this is observed in the evolution of the SMF, 
where, at high redshift the model predicts a later assembly of 
galaxies, exacerbating the underproduction of massive galaxies 
at higher redshifts. 
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Fig. 6. Comparison between GECO's quasar luminos- 
ity functions at different redsh i fts a nd observations 
by iHopkins. Richards. & Hernquisll |2P07. The magenta 
shaded region indicates the quasar luminosity function of 
iHopkins. Richards. & Hernquistll20()7l while GECO predictions 
are represented with a solid black line. 

It is worth to note, however, that such a comparison must be 
taken with caution before deriving firm conclusions. The errors 
in the stellar mass estimates can strongly affect (in the sense of 
increasing) the counts in the bright-end where the mass function 
is very steep. An improved approach would be to convolve the 
model stellar mass function with the error distribution in t he stel- 
lar mass given by the observatio ns ( Kitzbichler & Whitq, l2006t 
IStringer et alil2009tlCattaneo et aL,2008.) . Of course the uncer- 
tainties in the stellar mass estimates increase with redshift and 
this can influence the amount of the observed evolution of the 
stellar mass function. 

Finally, it is not unlikely that the derivation of the stellar 
mass function at the faint-end might have been seriously influ- 
enced by incompleteness effects in the source selection and have 
produced a spurious evolution at low masses. More extensive 
deep imaging by Spitzer will be carried out on large areas dur- 
ing the warm mission phase, which will offer an opportunity to 
address these issues shortly. 

6.2. Quasar luminosity function 

The gas accretion onto black holes described in Sect. 3.5.3 trig- 
gers the AGN activity whose bolometric luminosity associated is 
given by Lbh (eq.|2T]l. Although the accretion rate in our model 
is not limited by the Eddington rate, we limit the quasar lumi- 
nosity to lie below this limit: 

Lqso = ma\{LBH, LEdd) ■ (31) 

In order to test the consistency of the feedback from AGN ac- 
tivity on the galaxy population we compare the quasar luminos- 
ity function predicted by GECO with observational estimates. To 
avoid to introduce further uncertainties due the obscuration fac- 



tor in the optical bands, we compare our predictions with bolo- 
metric luminosity functions. We refer to the compilation from 
Hopkins, Richards, & Hernquist (2007), obtained by collecting 
a large multi-wavelength dataset, from the optical, soft and hard 
X-ray, near- and mid-IR. The authors give obscuration-corrected 
luminosities, aimed to represent the intrinsic quasar population. 
In Figure|6]we show our predicted quasar luminosity function up 
to z=2. The model gives an acceptable description of the evolv- 
ing quasar population at low and intermediate redshift, but at 
higher redshift (panel at z=2) it underpredicts the number of 
luminous quasars. Varying the efficiency of the accretion onto 
black hole {Jbh and kACN) does not solve the problem, because 
the match at high redshift would be obtained at the expense of 
the agreement of the local black-hole mass function. This may 
be an indication that the parametrisation of the AGN activity 
used in the pr esent work is not co mpletely adequate, as already 
suggested by iMarulli et al.l (12008 ), and an accretion efficiency 
increasing with redshift may be requh'e d (as sugges ted even by 
recent hydrodynamical simulations, see I Viola et al.l l2008). to fit 
simultaneously the low and the high redshift luminosity func- 
tion (Marufli et al., 2008; Lapi et al., 2006). We also suggest that 
higher accretion efficiencies at high redshift may be reached also 
with an enhanced cooling efficiency at early times. This alterna- 
tive may also provide a way to form galaxies more efficiently 
at high redshift, hence, improving the match of the stellar mass 
function evolution. 

6.3. Formation Histories of tiie Stellar Populations 

In order to get a deeper insight into the process of the stel- 
lar build-up in our modelled galaxies, we investigate here their 
detailed star formation histories (SFH). We derive the SFH a- 
posteriori for each present-day galaxies, taking into account its 
baryonic tree and summing up all the stellar populations formed 
inside its progenitors. We compute the SFHs over 52 time-steps, 
our pre-defined grid of times used for the merger tree. We show 
in Figure |7] as an example, the SFHs for 4 different realizations 
of the central galaxy of a halo with mass 5 x lO'^M©, i.e. very 
similar to the MiUcy-Way halo. We indicate the stars formed in 
the quiescent mode with a blue shaded histogram and with the 
red one the stars formed during the starburst phase. The black 
envelope is the sum of both contributions. In each panels, the 
host halo mass and stellar mass are indicated, together with the 
parameter r, that is the ratio between the bulge and the total stel- 
lar mass. 

In the case of such a Milky- Way-like halo, the majority of 
the central galaxies are spirals (~ 60%). However, the mix of 
galaxy properties turned out to be rather heterogeneous. In the 
lower-right panel we see an example of quite smooth star for- 
mation history with a steep rise of SFR at early times followed 
by a smooth decrease. In the two upper panels we see that more 
than one peak of quiescent star formation are present, indicating 
that star formation has occurred in different progenitors of the 
final galaxy. In other cases, the starburst can occur at more re- 
cent times (z ^ 0.5 in the lower-left panel) and is more efficient 
in consuming all the gas, hence the galaxy evolves passively un- 
til present. In all the cases represented starburst events are very 
frequent at early times, even in galaxies that at the present time 
have a late-type morphology. 

In order to give an indication of the age of a galaxy, we define 
its formation redshift as the redshift when the galaxy has formed 
half of its present-day stellar mass. We indicate it in Fig. [T] as 
a vertical dashed line. Galaxies with a smooth, quiescent star 
formation tend to form their stars later, while galaxies having 
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Fig. 7. SFHs of stellar populations for different realizations of the central galaxy of a Milky -Way like halo, having masses equal to 
5 X IO'^Mq. These SFHs sum up the contributions to the final galaxy stellar mass by all the stars present in the galaxy at z=0. Blue 
shaded histograms show the star formation occurring in the quiescent mode, the red ones represent the bursty mode, while the black 
envelope is the total SFR. The vertical dashed line indicates the formation redshift of the galaxy. The final stellar masses in each 
realizations are indicated. 



experienced a more bursty star formation history, thus having an 
elliptical morphology, tend to have higher formation redshifts. 

In Figure [8] we show the SFH of the central galaxy of ha- 
los with different mass, averaged over 100 realizations of the 
same parent halo tree. Halo masses range from ~ 10"Mq up to 
~ 10''*Mq. Note that, while in Fig. |7] the SFHs refer to single 
realizations of one halo, each panel here refers to the average 
over 100 realizations of the same halo, hence SFHs result to be 
smoother Two trends with the halo mass can be noticed. The first 
is that, passing from low-mass to high-mass halos, the starbursts 
become the dominant mechanism of star formation. The second 
is that the star formation rate for galaxies living in high-mass 
halos peaks at higher redshift compared to that in lower mass 
objects. Therefore, the associated formation redshift for massive 
objects is higher. In fact, the increasing importance of the star- 
burst mode, together with the AGN feedback connected to it, 
lead to a fast consumption of gas at early times, by preventing 
the gas to cool and inhibiting further star formation. The con- 
sequence is a lower rate of star formation at later times, hence 
higher average formation redshifts for galaxies living in high- 
mass halos. 

Since high-mass galaxies live in high-mass halos, these re- 
sults indicate that stars in massive galaxies are on average older 



than in their less massive counterparts. Nevertheless even in 
massive galaxies the star formation is not completely quenched, 
but a long tail of star formation at low redshift still occurs. 
This behaviour of our SFHs appears to be in general agree- 
ment with the evidence for a downsizing pattern of galaxy 
evolu t ion, an effect discussed since long time (ICowie et all 
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19961: Gavazzi & Scodeggiol 119961 : iFranceschini et al.l 

Thomas et al. I2005h . We will discuss further the physical rea- 
sons for this in Sect. ^ 

6.4. The Cosmic Star Formation Rate Density 

The evolution of the comoving star formation rate density for 
the whole galaxy population is shown in Figure |9] (left panel). 
The model prediction is shown by the black line and is com- 
pared with data from the compilation by Hopkins (2004), for 
UV, optical and radio tracers, and from lRodighiero et al.i (i2009l) 
for the IR tracers. Different colours encode different tracers for 
the SFR: cyan for star formation derived from the UV contin- 
uum, green for optical tracers (O ii. Ha and HyS), magenta for 
data using radio tracers and red for SFR estimates derived from 
the IR luminosity. The latter are believed to offer particularly 
reliable and robust estimates of the star-formation rate because 
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Fig. 8. Averaged SFHs for central galaxies of 4 different halos, with masses ranging from ~ 10"Mo to ~ IO'^^Mq. The average are 
computed over 100 realizations. As in figure|7] blue shaded histograms show the star formation occurring in quiescent mode, the 
red ones represent the bursty mode, while the black envelope is the total SFR. The host halo mass and the averaged final stellar mass 
are indicated. 



they trace rather directly the emission by you ng stellar popula- 
tions and are free of dust extinction effects jKennicuttl Il998t 
iFranceschini et al.L 120011: lElbaz et al.L l2002h . Th e blue star at 
z^7 is the determination of the SFR density of iBunker et alj 
(12009') through HST-WFC3 data with the z' -drops technique, 
effective over the redshift range 6.7<z<8.8, as indicated by the 
X-axis bar 

The model provides a remarkably good fit to the data over 
most of the redshift range. It should be stressed, in particu- 
lar, the nice agreement of the model with IR data, showing a 
fast rise from z = to z = 1 and then a wide plateau up to 
z ~ 4. Locally (z < 0.2) there is an indication that our model 
produces too many stars. It appears consistent with UV deter- 
minatio ns, but it keeps h igher with respect to the determina- 
tion of iRodighiero et all (l2009h using the Spitzer 24fim lumi- 
nosity function. Even at very high redshift (z^ 7), the model 
is marginally consistent with the observations, especially if we 
take into account that the HST-WFC3 data can be considered as 
a lower limit due to incompleteness effects. We note also that 
the WMAP5 model predicts a much lower SFR at high redshift, 
(z > 2) leading to a better agreement with the z^ 7 data. Contrary 
to what found in Somerville et al. (2008), with the new cosmol- 
ogy we find a better agreement with data in the SFR density, but 
a worst one in the stellar mass assembly, as noted above. 



The high level of star formation rate at redshifts z ^ 3 - 5 
found in the model was already mentioned in the previous Sect. 
On one hand, the quiescent mode of star formation is very effec- 
tive, due to the high cooling efficiency in denser environment. 
Indeed, the peak in the quiescent mode of star formation (blue 
shaded histograms in Figs. [8]l is always at z ~ 2. On the other 
hand at high redshift the star formation mainly occurs in star- 
bursts, triggered by the high merging rate, which enhances the 
SFR at early cosmic times. The exhaustion of cold gas available 
explains the sharp later decline at lower redshifts. 

In the right panel of Figure |9]the contribution to the total star 
formation rate density (black solid line) is splitted according to 
the mass of the hosting halo: the short-dashed red line shows the 
contribution from galaxies living in high-mass halos (M > 3 x 
10 '^Mo), while the green long-dashed line shows the one from 
galaxies living in less massive halos. The time dependence of 
the star formation rate in the two cases has a similar behaviour: 
a slow increase at high redshifts, a phase of peak activity, and 
finally a quick decline to the present time. However, in addition 
to this general pattern, we notice a systematic shift in the SFR 
history between the two halo populations, with the high-mass 
halo activity being shifted at higher redshifts, z ~ 2 - 4, while 
it is more concentrated around z = 1 - 2 in the less massive 
population. 
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Fig. 9. Left panel: SFR density as a function of redshift. Solid black line shows the GECO predictions. The shaded area and coloured 
points are the SFR obtained from observations at different wavelengths. The cyan shaded region shows the SFR derived from UV 
tracer s, green region means that observations are in the optical, and magenta dots in the radio (from the compilation of Hopkins| 
red squares refer to infrared de terminations, from Rodighiero et al. (2009.) and the blue star is the z -1 determination from 
HST-WFC3 data dBunker et al.L l2009h . Right panel: SFR density spUtted in the contribution by objects residing in different halo 
masses. Solid black line shows the total SFR density, red short-dashed shows the contribution from galaxies living in high-mass 
halos {Mhaio > 3 X IO'^Mq), while the green long-dashed curve represents the contribution from galaxies in less massive halos 

(MHalo < 3 X lO'^Mo). 



We then confirm our previously found results: galaxies re- 
siding in high-mass halos form their stars at earlier epochs, and 
then their star formation rapidly slows-down, while galaxies in 
less massive systems form stars at an enhanced rate even at re- 
cent times (z < 1). 

6.5. Star Formation and Mass Assembly 

We have discussed in the previous Sect, an important aspect of 
our model, its prediction that massive systems form their stars 
at earlier times then less massive objects. This however con- 
cerns the birth times of the galaxy's stellar populations, while 
the times of mass assembly may be very different when stars are 
formed at high redshift in a number of distinct progenitors and 
are assembled at more recent times. 

In Figure[TO]we compare the star formation history with the 
"assembly history" for single realizations of various parent ha- 
los. In the upper boxes of each panel we show the SFH of all stel- 
lar populations contributing to the final galaxy, in analogy with 
those reported in Fig.|7] The SFH is expressed here in terms of 
the fraction of the final stellar mass formed in each time-step (in- 
stead of the SFR), for an easier quantitative comparison with the 
assembly history. In the bottom boxes we show the mass assem- 
bled into the main progenitor at each time-step, where the main 
progenitor identifies at each redshift the progenitor of the galaxy 
which survives until the present day. Here, red histograms indi- 
cate the amount of stars added to the bulge in a given time-step, 
while blue histograms show the stars added to the disc. Although 
in a few cases the history of stellar assembly is very similar to 
that of star formation, indicating that star formation took place 



mainly in one single object, in general the two paths may be sig- 
nificantly different. 

In most of the cases the assembly of bulges, due to merg- 
ers, is related to starbursts, although in some systems, mainly in 
merging episodes at low redshifts, the mass assembly is not as- 
sociated to any event of star formation because the gas content 
is already consumed (dry merging). Substantial dissipationless 
merging may lead to very different formation times for the pro- 
cess of star formation and mass build-up. Although dry mergers 
in the local universe are observed, their frequency and their role 
in galaxy evolution are still poo rly constrained by observations 
(Ivan Dokkumi l2005t iBeU" et al.L I2OO61) . Anyway galaxy forma- 
tion models predict that th ey are a crucial proce ss in the forma- 
tion of local ellipticals ("see lCattaneo et al]l2010h . 

We define the redshift of assembly as that corresponding to 
when half of the final stellar mass is assembled in one single ob- 
ject, and we show it in Fig. [10] By definition, the assembly red- 
shift is always smaller than the formation one, they equal each 
other when almost all the stars are formed in the main progenitor 

We report in Figure [TT] the dependencies of the mean star 
formation redshift (upper panels) and the redshift of mass as- 
sembly (lower panels) as a function of the hosting halo mass 
(left-hand panels) and of the galactic stellar mass (right-hand 
panels). Error-bars indicate the standard error of the mean over 
100 MC realizations of the same host halo. Note that in order to 
derive the redshift of formation (or assembly) as a function of 
stellar mass we weight each galaxy mass by the number density 
of its host halo, because any stellar mass may belong to a certain 
range of halo masses. 
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Fig. 10. Comparison between stellar mass formed (upper panels) and assembled into the main progenitor (lower panels) for different 
model galaxies, living in different halos. In the upper panel red and blue histograms represent, respectively, bursty and quiescent SF, 
while in the assembly panels red and blue stay for the mass assembly of bulges and discs respectively. 



The mean formation redshift shows a clear increase with 
both halo and stellar mass. This effect is more pronounced for the 
central galaxies than for the whole population including satel- 
lites. We have already shown in the previous Sects, how galaxies 
living in high-mass halos formed their stars earlier and at a faster 
rate than in low-mass halos. Here, we demonstrate that the trend 
is true even with the galactic stellar mass, although it is weaker 
Note that the most massive halos (M/,a/o > lO'^^M©) host galax- 
ies with the higher formation redshift (z/o™ ~ 2.5, top left panel 
of the figure), but since these halos are very rare, the mean for- 
mation redshift of the more massive galaxies {Msiar ~ lO'^M©), 



which belong to a certain range of halo masses, keeps somewhat 
lower (~ 2, top right panel). 

The behaviour of the redshift of mass assembly is more dif- 
ficult to interpret. We found that, on average, galaxies living in 
high mass halos assemble their stars earlier Indeed, the merger 
time-scales due to both dynamical friction and random collisions 
among satellites increase with the halo mass. We already men- 
tioned in ^3.6.11 that the dynamical friction time increases with 
the ratio between the mass of the halo and that of the satellite 
(eq. |23] ). Moreover, since the merger probability is strongly re- 
duced in high-velocity encounters, the merger rate for satellites 
decreases with the ratio Vc/vgai (eq.l24l). and becomes very low 
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Fig. 11. Mean formation redshift (upper panels) and assembly redshift (lower panels) of galaxies as a function of their host halo 
mass (left-hand panels) and stellar mass (right-hand panels). Black squares indicate the global population of galaxies, including 
central and satellites, while red open circles stay for central galaxies. Error-bars indicate the standard error on the mean over 100 
MC realizations of the same halo. 



for massive cluster-size halos. Nevertheless, when the redshift 
of mass assembly is expressed in terms of the galactic stellar 
mass, the downsizing behaviour is no more recovered: the as- 
sembly time keeps approximately constant around Zass ~ 1 for 
low and intermediate masses and decreases in the highest mass 
end, with Zass ~ 0.7 for Mstar ~ lO'^M©. These stellar mass val- 
ues correspond to those of the CD galaxies in local rich clusters 
or groups. The combined effects of the formation redshift and 
the assembly redshift suggest that the only way to form the most 
massive galaxies in the local universe is in the form of dry merg- 
ers at recent times. Indeed, massive galaxies cease to form stars 
quite early and, without such events, they are not able to grow 
their mass above Mstar - IO^Mq. 



In conclusion, our model expectation is that, although merg- 
ers and collisions are rare in massive halos, they do happen, and 
their effect is crucial in the formation of the today most mas- 
sive galaxies. Refined observational estimates of the galaxy mass 
functions at the high mass-end all over the z=0 to 4 redshift in- 
terval (on top of those already apparent in Fig. |5]) will clarify 
how far this expectation of a substantial dry merging process is 
born out by the data. 



7. Discussion and Conclusions 

We have presented a new semi-analytical model of galaxy for- 
mation, the Galaxy Evolution Code, GECO which appears to 
reproduce several key statistical properties of local and high- 
redshift galaxies. 

GECO is based on a state-of-the-art Monte Carlo algorithm 
for the representation of the dark matter halo merging history, 
based on the Extended Press-Schechter formalism. GECO in- 
cludes detailed implementations for gas cooling, star formation, 
feedback from SN and galaxy mergers, due to both dynamical 
friction and random collisions. Moreover, the parallel growth of 
BHs is followed in time and the subsequent AGN feedback is 
modelled. 

We specifically tested our results directly on the observables 
involving the stellar mass and star-formation rate more then lu- 
minosities, as usually done by other published models. This is 
motivated on one hand by the fact that the stellar mass functions 
are the most direct outcome of the model. On the other hand, 
stellar masses in galaxies have recently become a rather straight- 
forward observable thanks to rest-frame near-infrared data by the 
Spitzer Space Telescope, directly probing the stellar mass con- 
tent in high-redshift galaxies. At the same time, Spitzer is also 
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probing with deep far-infrared photometric imaging the rate of 
stellar formation in distant objects. Therefore, we can compare 
the outcomes of our model with the most reliable set of observ- 
ables and a minimal number of free parameters. We thus believe 
that the basic physics of the ACDM hierarchical clustering con- 
cept of galaxy formation can be tested by us in a very effective 
way. 

The main results obtained in the present work are sum- 
marised in the following. 

1 . The local stellar mass function result s in a remarkably good 
agreemen t with t he determination of ICole et al.l (1200 lb and 
iBell et all (l2003h (Fig. At the high-mass end the total 
mass function is dominated by the contribution of bulges, 
while discs dominate at the faint-end. When the total stel- 
lar mass is splitted into the contributions of early-type and 
late-type galaxies, the former populate the bright-end side, 
while the latter mainly contribute at low masses. The num- 
ber densities of the two morphological types cross each other 
at Msiar ~ 3 X IO^'Mq/;"', as observed. Although we repro- 
duce the general trend of the morphological mass functions, 
our model fails in matching the low-mass end of spheroids, 
showing an excess of low-mass systems. Likely, this is due to 
an oversimplification of the satellite population, since satel- 
lite galaxies loose their hot gas reservoir as soon as they are 
incorporated in a more massive halo and the star formation 
is quenched soon after. 

2. The co-evolution of galaxies and BHs is modelled following 
the prescriptions of Croton et al. (2006). A first mode of ac- 
cretion onto BHs considered is the so called 'radio-mode', 
that inhibits the quiescent star formation, while the second 
one is the 'QSO-mode', that is triggered only during galaxy 
mergers and constitutes a major channel of BH accretion. 
As a consequence of mergers, a starburst is induced as well, 
feeding the galactic bulge component (and destroying the 
disc in the case of a major merger). This leads to a parallel 
growth of BH and bulges with the two masse s very well cor- 
related, in agreement with observational data ('H aring & RixL 
|2004) and to the local black-hole mass function in a re- 
mark able agreement with the observations of (IShankar et al.L 
12004 . 

3. We compare the stellar mass functions resulting from the 
model with various observational determinations up to z ~ 
3.5 (Fig. |5]) and found a reasonably fair agreement over the 
whole redshift range considered. Nevertheless, the observed 
ratio between the evolution of the faint- and the bright-end 
of the stellar mass function is not very accurately repro- 
duced: there is too much evolution in the model at the bright- 
end and too little at the faint-end compared to observations. 
Various sets of observables indicate a large increase in the 
number density of low-mass objects between z ~ 2 and the 
present-day and a lower rate of evolution for massive ob- 
jects. However, we mentioned that the completeness and ro- 
bustness of the observational mass function are to be proven 
there, before claiming more definite conclusions. In the case 
of WMAP5 cosmology we observe a delay in the formation 
of cosmic structure, that lead to a further reduction of high- 
mass systems. 

4. The bolometric quasar luminosity functi on is compared with 
iHopkins. Richards. & Hernquisj (l2007h . showing a good 
level of agreement at low and intermediate redshift, but a ten- 
dency to underpredict the number of bright quasars at high 
redshift. A mechanism for enhancing the cooling rate at high 
redshift might simultaneously increase the fuel for BHs and 



enhance the SF at high redshift, as seems to be required to 
improve the match with the stellar mass functions. 

5. The integrated star formation rate density (Fig. |9j shows 
an high level of star formation at high redshifts, a peak at 
z ~ 1.5-3 and then a sharp decline below z ~ I- When 
compared with the determination of the SFH derived using 
various tracers (UV, optical, radio, IR), our predictions are in 
very close agreement with these observations. At very high 
redshift (z - 7) our model is able to correctly reproduce the 
recent determination of the SFR density by HST-WFC3 data. 

6. We analysed in detail the SFH in simulated sets of galax- 
ies, to gain insights into how the model treats star formation 
and how it depends on the galaxy or halo mass. We com- 
puted the averaged SFH for the central galaxies living in ha- 
los of different sizes (Fig.[8]l. We identified two main trends 
with halo mass. First, going to high-mass systems, the con- 
tribution of the starburst mode to the total SF becomes in- 
creasingly important, and indeed predominant in very high 
mass objects. Second, the formation redshift, defined as that 
when half of the present-day stellar mass is formed, in- 
creases, leading to older stellar populations in massive sys- 
tems. Hence galaxies in our model form their stars follow- 
ing a downsizing pattern, consistent, for inst ance, with the 
dating of stellar populations in local galaxies dThomas et al.L 
l2005h . The naive expectation of early versions of hierarchi- 
cal galaxy formation models was that, since massive halos 
are assembled later than their lower-mass counterparts, the 
most massive galaxies, hosted in the largest halos, should 
form their stellar content at the same late cosmic time. 
Actually, as shown in Fig. (TT] this is clearly not the case 
in our refined model. According to it, downsizing in star- 
formation is an intrinsic feature of semi-analytical mod- 
els (see also iNeistein. van den Bosch. & Dekelll2006l) . The 
present-day massive galaxies were formed through the as- 
sembly of a number of smaller progenitors, that collapsed at 
high redshifts from the highest density peaks of the primor- 
dial density field. Acc ording to this scenari o, also named bi- 
ased galaxy formation dPekel & SiUd,ll986l) . bright and mas- 
sive systems started to form stars early on. This feature is a 
natural outcome of the merger tree formalism: progenitors 
of high-mass halos fall below the resolution mass imposed 
to the merger tree after several time-steps back in cosmic 
time, so the leaves of the tree are found at high redshift. On 
the contrary, smaller systems, closer to the resolution mass, 
take only a few time-steps back to reach this minimum mass. 
Note that this not merely a computational artifact, but indeed 
corresponds to the fact that we expect quite negligible star 
formation to have occurred below such threshold mass, on 
consideration of the SF quenching by the UV photoionizing 
background (Sect. I3.5.2l i. Since baryons are put into halos 
starting from the leaves, in high mass halos star formation 
took place at early times. Moreover, at high redshift both 
mechanisms of star formation were more effective. Thanks 
to the efficient cooling of the gas, the quiescent mode of star 
formation occurs at enhanced rate. Also, the frequency of 
mergers at these times is high, allowing an efficient conver- 
sion of gas in stars via a starburst. The exhaustion of cold gas 
is, therefore, very rapid, leaving the galaxy devoid of fuel, 
and preventing further star formation to occur 

A further reason for the star formation quenching in mas- 
sive systems may be ascribed to AGN feedback. In order 
to check the importance of such mechanism, we run a test 
simulation with the AGN emission switched off. The com- 
parison between the two versions of the model are presented 
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in Appendix A. We found that AGN feedback has some ef- 
fect in increasing the average age of stellar populations in the 
most massive galaxies, as it was poin t out by previous works 
dCroton et all l2006t iDe Lucia et all 120061 ICattaneo et all 
I2008h . However, in our model the downsizing trend is still 
obtained, although slightly weaker, even in the absence of an 
AGN. 

Therefore, in GECO, the AGN feedback has not a dramatic 
effect in producin g th e local galaxy properties (see also 
iMenci et alJ l2006l and iMonaco. Fontanot. & Taffonil l2007h 
and it is not the main reason for the downsizing pattern of 
galaxy evolution. In our modelling, the latter is rather due to 
a higher efficiency of SF back in time for the most massive 
galaxies in rich environments, as explained above. Although 
we have shown that the downsizing is quite a natural feature 
of our model, it is not the case for all the published models 
and some of the m fail in repr oducing such trend with mass, 
as it is shown in lFontanot et al. (.2009) . 
7. Finally, we compared the star formation history of simulated 
galaxies with the mass assembly history, that is the history 
of the mass assembled into the main progenitor at each time- 
steps (Fig. [Toll. As expected, the two processes may occur on 
very different time-scales, especially in high-mass systems, 
where the star formation took place at high redshift in many 
distinct progenitors, which assemble at low redshift. The late 
assembly of these systems occurs in the majority of the cases 
through dry mergers, i.e. mergers between spheroidal sys- 
tems with little or no gas. In these cases, the merger is not ac- 
companied by any event of star-formation, neither quiescent 
nor bursting. An intriguing question is then whether galax- 
ies in the model assemble their stellar mass in a downsizing 
way, as it occurs for the star formation. Possible downsizing 
effects in the mass assembly wer e indicated by several au- 
thors (Bundv et al., 2006; Cimatti . Daddi. & RenziniL l2006t 
iHopkins et al.1, 12007 ) as inferred from the lack of evolution 
in the high-mass end of the mass function and by the evolu- 
tion of the blue-to-red galaxy crossover mass (the mass for 
which the early-type mass function intersects that of late- 
types). In our model, we found that although the assembly 
time shows a shallow dependence on the host halo mass, 
on average galaxies living in massive halos assemble their 
stars before galaxies in less massive hosts (Fig. fTTT i. Instead, 
the assembly times are almost constant with stellar mass, 
and decrease for very-high mass systems {Msmr > 1O"M0), 
hence leading to an upsizing trend with time in the high- 
mass end. In our modelling, dry mergers are the main reason 
for this late assembly of massive galaxies. These findings 
agree with other previously published semi-analytical mod- 
els (De Lucia et al., 2006; Cattaneo et al., 2008, 2010). The 
importance of dry mergers in the formation of the most mas- 
sive galaxies that we observe in the nearby universe, is im- 
plied also by the re cent observations of the size evolution o f 
massive spheroids (Truiillo e t all 120071; ICimatti et alll2008l) . 
Indeed, if such compact gal axies are absent i n the local uni- 
verse as suggested by Tr uiillo et al.l (l2009h . the expected 
mechanisms that move the high-redshift compact galaxies 
to the local relation are dissipationle ss mergers. Anyway, 
other works (Vale ntinuzzi et all 1201 Oh suggest that such su- 
perdense galaxies in the local universe are not as rare as pre- 
viously claimed. 

How far this might be incongruous with the observational 
indication of downsizing in mass assembly mentioned above 
will become clear with furth er observational confrontation. 
For example, ICattaneo et al.l (i2008) argue that the upsizing 



in mass assembly can coexist with a downwards trend in the 
transition mass, which then turns out to be a poor indicator 
of downsizing. 

In conclusion, GECO presents an encouraging level of agree- 
ment with a wide range of observational data, at low and high- 
redshifts. We focus on comparing GECO with data on the two 
main phases of the galaxy formation process, that are the star 
formation and the mass assembly. On one hand, we confirmed 
that the observed downsizing in star formation is natural part 
of our scheme of hierarchical growth of structures. On the other 
hand, the stellar mass assembly process remains more difficult to 
understand from both a theoretical and an observational point of 
view. The times of galaxy assembly in our model, related with 
both the galaxy merger time-scales and the star formation ef- 
ficiency, strongly depend on the details of the implementations 
of galaxy dynamics (dynamical friction, satellite collisions and 
tidal stripping): further work in this sense remains to be done 
in order to have a deeper insight into the galaxy assembly pro- 
cess, and certainly a comparison with N-body simulations will 
be helpful. 

The most striking conclusion is that, despite the simplicity of 
the prescriptions adopted, and the small number of free param- 
eters used, the main features of the evolving galaxy population 
are reproduced. In particular, the AGN feedback is needed only 
to improve the match of the local stellar mass function, but its 
effect on stellar ages is not determinant. 

We have to keep in mind in any case that this paper includes 
just a preliminary and partial confrontation with the data, show- 
ing at least no obvious clash. Much more extensive analyses and 
tighter constraints will be obtained as soon as refined data on the 
evolutionary mass functions and stellar birthrates will become 
available. 
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Appendix A: Effect of AGN feedback 

Several recent works (ICroton et al.L 120061; iMalbon et al.L l2007t 
ISiiacki et all \200l h advocated the AGN feedback as the only 
mechanism able to quench star formation in order to match the 
local abundance of massive galaxies as well as their old ages and 
red colours. Here we check the relevance of such feedback in 
reproducing the observations in GECO, comparing our fiducial 
model with a version in which the AGN feedback (both in quasar 
and radio modes) is switched off. In Figure IATTI we can see the 
effect of the AGN in the stellar mass function. As expected, the 
presence of the AGN influences only the bright-end, reducing it 
by ^ 0.15 dex, both in the local universe (upper panel) and at 
high redshift (lower panel), hence leading to the same amount 
of evolution in the stellar mass function as in the fiducial case. 
We note that given the overabundance of massive galaxies by the 
same amount at all the redshift, in the model without AGN the 
comparison with the observational mass function at high redshift 
results in a closer agreement with respect to the fiducial model. 
Figure IA.2I shows that, obviously, in the model without AGN 
the star formation keeps lower at all redshifts, with an indication 
that at very high redshift (z > 4) the difference between the two 
models tends to disappear, meaning that the importance of AGN 
feedback increases at late cosmic times. 

A further issue to be handled is the relevance of the AGN 
in producing the downsizing trend in the stellar ages. In Figure 
IA.3I the trend of the formation redshift with mass is shown for 
both models. Although the AGN is effective in increasing stel- 
lar ages, the n et effect is act ually very modest. At variance with 
other models dCroton et al.l 12006. for instance) the downsizing 
trend is still recovered, leading to the conclusion that the pres- 
ence of AGN is not the main reason for the onset of the down- 
sizing pattern for galaxy evolution in our model. 

Therefore, in our model the AGN has the effect of reducing 
the star formation activity in the high-mass objects, but the 
amount of reduction is perhaps less striking than previously 
claimed (Croton et al., 2006: Malbon et all 120071: iMarulh et al.[ 
[2008: Lagos, Cora, & Padilla, 2008). Indeed, the values for the 
efficiencies of black hole accretion adopted here, are somehow 
smaller than those adopted by the models mentioned above, but 
we found that they are suitable to obtain the very good match of 
the local relations. 
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Fig.A.l. Upper panel: comparison between the stellar mass 
function at z=0 of our fiducial model (dashed black line) with the 
model without AGN feedback (solid purple line). Lower panel: 
the same as the upper panel but for z=2.5. Observational data at 
the given redshift are also indicated by the shaded region. 
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Fig. A.2. Comparison between the SFR density as a function of 
redshift in the model without AGN feedback (soHd purple line) 
and the fiducial case (dashed black line). 
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Fig. A.3. Comparison between the formation redshifts as a func- 
tion of stellar mass in the model without AGN feedback (solid 
purple Une) and the fiducial case (dashed black line). 
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Appendix B: Effect of cosmic reionization 



As explained in Sect. 13.5.21 our treatment for the cosmic reion- 
ization is based on the resuhs previously found byl Benson et al.l 
(l2003h . Anyway, it has been claimed, on the basis of recent hy- 
drodynamical simulations, that the effect of the photoionizing 
background on st ructure formation is indeed weaker than pre- 
viou sly assumed dHoeft et al.L l2006t lOkamoto. Gao. & Theunsl 
12008 ). Taking into account such results, we test the effect of our 
assumption on the limit circular velocity for the reionization. In 
Figure IB. II we show a comparison between our fiducial model 
with a model in which the circular velocity limit for baryon cool- 
ing has been lowered down to 25 km/s. To allow the formation 
of low-mass halos, we used for this model a merger tree with 
a resolution of IO'^Mq. We see that a lower limit on the circu- 
lar velocity has the effect of slightly increasing the number of 
low mass galaxies. At the same time, the number of massive ob- 
jects is reduced, since a greater amount of gas has been locked 
within satellites and is not available for star formation inside the 
brightest systems. In any case, the effect is rather moderate and 
the prediction of the fiducial model could be accommodated, for 
instance, by increasing the effect of SN. 

Therefore, we conclude that although we might have over- 
estimated the effect of reionization in our fiducial model, the 
main results found in the previous sections keep unchanged if 
we allow a milder dependence of galaxy formation on the cos- 
mic reionization. 
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Fig. B.l. Upper panel: comparison between the stellar mass 
function at z=0 of our fiducial model (dashed black line) with 
a model with a reduced effect of the cosmic reionization (solid 
yellow line). Lower panel: the same as the upper panel but for 
z=2.5. Observational data at the given are also indicated by the 
shaded region. 



